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ABSTRACT 

Most observed extrasolar planets have masses similar to, but orbits very different from, 
the gas giants of our solar system. Many are much closer to their parent stars than 
would have been expected and their orbits are often rather eccentric. We show that 
some of these planets might have formed in systems much like our solar system, i.e. 
in systems where the gas giants were originally on orbits with a semi-major axis of 
several au, but where the masses of the gas giants were all rather similar. If such a 
system is perturbed by another star, strong planet-planet interactions follow, causing 
the ejection of several planets while leaving those remaining on much tighter and 
more eccentric orbits. The eccentricity distribution of these perturbed systems is very 
similar to that of the observed extrasolar planets with semi-major axis between 1 and 
6 au. 

Key words: Celestial mechanics, Stellar dynamics; Binaries: general, Planetary sys- 
tems 



1 INTRODUCTION 



Since their fir st discovery (jWolszczan fc Fraill Il992l; 

iMavor fc Quelozl Il995l ). more than 300 extrasolar planets 
have been found. In Fig. 1 we plot the semi-major axes and 
eccentricities of the planets de tected using the radial veloc- 
ity method as of October , 2008 (|Butler et al.l200rj ; [Schneiderl 
2008; Tamuz et al.|[2008l ). As can be seen from the figure, the 
spread in eccentricity and separation is very large. A major- 
ity of the detected extrasolar planets have masses similar 
too, or larger than, those of the gas giants in our solar sys- 
tem. Most planet formation models predict that such mas- 
sive planets can only form outside the so-c alled snow line 
(situ ated at 3 au around a solar-mass star ( Ko kubo fc 
1998)) but as can be seen from Fig. 1, most of the detected 
exoplanets have orbits tighter than this value. Thus, the or- 
bits of most of the observed exoplanets must have shrunk 
considerably since their formation. The most efficient mech- 
anism behind t his is widely believed to be disk migration 
|Lin et al.lll996l ). However, disk migration almost exclusively 
produce planets on circular orbits, which does not agree with 
observations. Thus, in order to explain the observed eccen- 
tricities an additional mechanism must be at work. 

The most popular such mechanism is scatterings due 
to strong planet-planet interactions which can explain the 
observed eccentricities. Several different models which re- 
produce the observed eccentricities very well have been sug- 
gested and most likely the observed sample is created from 
a combination of these and perhaps others. It may, for 
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example, be that planets come very close to each other 
while undergoing migration in the disk (|Moorhead fc Adamsl 
l2005h or that the planetary orbits are kept stable by the 
disk and as it evapo rates the system becomes unstable 
|Thommes et al.l 120081 ) . Another possibility is that many 
planetary systems are initially too tightly packed, lead- 
ing to that they become unstable on a time-scale of mil- 
lions to several hundred million years and undergo strong 
planet-planet in teractions (M arzari & Wcidcnschilling 2002; 
Barnes fc QuinrJ |2004|; [Ford et al.l 120051 ; IChatteriee et al] 
20071 ; I Juric fc Tremainel(2007l ). 



It is however also possible that strong planet-planet in- 
teractions are triggered in long-term stable planetary sys- 
tems (like, for example, our own solar system) by exter- 
nal perturbations. Two examples of such are nearby pass- 
ing s t ars in young stella r clusters ([Za kamska fc Tremaind 
120041 ; I Adams et all l200rj ; iMalmberg et all l2007bl ) and the 
effects of a stellar companion in a binary (see for example 
iHolman et al. 1 1 19971 ; iMazeh et al.lll997l : IHolman fc Wiegertl 
1999). In the latter case there are two distinctly different 
scenarios. If the planetary system formed in a primordial 
binary system, the orbits of the planets and the companion 
star are expected to be essentially co-planar. The evolution 
of s uch a system h a s bee n studied extensively in, for exam- 
ple, M 

arzari et al.l (|2005) . If the planetary system instead 
formed around an originally single star, which was later ex- 
changed into a binary in an encounter in a young stellar 
cluster, the orientation of the orbits of the planets is com- 
pletely random with respect to the orbit of the companion 
star. The evolution of the system is then very different from 
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Figure 1. Eccentricity, e, of the extrasolar planets detected using 
the radial velocity method plotted against their respective semi- 
major axes, a (in au). The crosses identify planets orbiting single 
stars, while the filled circles identify planets in binaries. The two 
open squares represent Jupiter and Saturn. The grey region in this 
plot encompasses the region in phase-space, which we consider in 
this letter. 



the co-planar case jTakeda fc Rasiol 120051 : iMalmberg et al.l 
l2007al ; iMarzari fc Barbierill2007h . 

In this letter we study how a companion star in a bi- 
nary, which formed through an exchange encounter in a 
young stellar cluster, can affect solar-system- like planetary 
systems. A solar-system-like planetary system is by us de- 
fined as a planetary system in which the gas giants are on 
long-term stable orbits wider than 5 au. We compare the 
eccentricity distribution of the resulting systems with the 
eccentricity distribution of the observed extrasolar planets 
at intermediate separations from their host stars. We find 
that they are rather similar if most of the planetary systems 
originally consisted of gas giants with rather similar masses. 



2 SAMPLE SELECTION 

We define planets at intermediate separations as those with 
semi-major axes in the range 1-6 au. The upper limit roughly 
coincides with the observational limit and is only a little 
larger than the semi-major axis of Jupiter in the solar sys- 
tem. Hence, 6 au roughly corresponds to where we would 
expect to start finding planets in unperturbed solar-system- 
like planetary systems. From our simulations of we find that 
it is very unlikely for planets to be scattered onto orbits 
tighter than about 1 au. This is not an exact limit, and we 
do see a few planets on even tighter orbits in our simulations. 
It is however clear that the planetary systems produced in 
our simulations will essentially only contribute to the ob- 
served extrasolar planets with semi-major axis greater than 
about 1 au and hence we only compare our results to the 
observed planets with a > 1 au. 



According to ICummind (|2004f ) the detection efficiency 
of planets using the radial velocity method decreases sharply 
with increasing eccentricity for planets with e > 0.6. This 
implies that the observed eccentricity distribution is not 
co mplete above an ecce ntricity of 0.6. However, according 
to lShen fc Turner! {2008) the decrease in detection effic iency 
is not as strong as that predicted by ICummind (|2004l ) and 
hence the observed eccentricity distribution is a good re- 
flection of the true eccentricity distribution. Nevertheless, 
in order to be certain that we avoid comparing our simula- 
tions to a biased sample we here only consider the planets 
in the shaded area of Fig. 1, hence those with semi-major 
axes between 1 and 6 au and with eccentricity less than 0.6. 
We have calculated the cumulative eccentricity distribution 
of this sample, and find it to be increasing approximately 
linearly between an eccentricity of and 0.6. 



3 PLANETS IN BINARIES 

In Fig. 1 we have identified the planets found in binaries 
(filled circles) from those found around single st ars (crosses) 
Pesidera fc Barbieril [20071 ; iTamuz et all 120081 '). It is clear 
from the figure that there is no obvious difference in the 
eccentricity distribution for planets in binaries with respect 
to that of planets orbiting single stars for the planets in the 
shaded region. It is however evident that the four most ec- 
centric systems are found to be in binaries. Since the number 
of systems is very low it is however too early to say if this 
is a real effect, or just a coincidence. 

Using a thorough statisti cal analysis of the o b served 
sample of extrasolar planets, iDesidera fc Barbieril (|2007T ) 
showed that there is no significant difference between the 
eccentricity distribution of planets in binaries and that of 
planets orbiting single stars below e = 0.6. Above e = 0.6 
there is possibly a slight excess of highly eccentric planets 
for the planets in (wide) binaries, although the statistical 
significance of this finding is, due to the low number of sys- 
tems, not rigorous. This is a very important result, since it 
leaves only two possibilities: 

(i) Planetary systems are not affected by the presence of 
a companion star in a binary, or 

(ii) the perturbation from the companion star in a bi- 
nary triggers the same mechanism as that which give rise to 
the observed distribution of eccentricities for planets around 
single stars. 

The separation of most observed planet-hosting bina- 
ries is between 100 and 1000 au and the companion mass is 
between 0.2 and 2.0 M© (|Desidera fc BarbierJl2007h . If most 
of the binaries with 100 < a < 1000 au are not primordial, 
but were instead formed in exchange encounters in young 
stellar clusters, the orientation of the orbital plane of the 
planets with respect to the orbital plane o f the companion 
star is completely random (see for example IMalmberg et al] 
|2007a|). In that case about 77 per cent of the systems will 
have an inclination between the orbits of the planets and the 
orbit of the companion star greater t han 39 . 2°. In such sys- 
tems the so-called Kozai Mechanism (|Kozaill 19621 ) operates, 
making the orbits of the planets more eccentric. This can, 
depending on the initial configuration of the system, trigger 



© 2008 RAS, MNRAS 000, \TM 



On the origin of eccentricities among extrasolar planets 3 



strong planet-planet interactions. For example, the four gi- 
ants in our solar system would, if put inside a binary with 
properties similar to those observed for planet-hosting bi- 
naries, undergo a phase of strong planet-planet interactions 
within a few million years, assuming that the inclination 
of the companion sta r with respect to the p lanets is large 
enough (see Fig. 2 in iMalmberg et al.ll2007bh . 

Strong planet-planet interactions leave almost no trace 
of the initial eccentricities of the planets and it is thus most 
likely the dominating process behind the observed eccen- 
tricities among planets in binaries. Hence, alternative (i) is 
not correct, leaving us with alternative (ii): the perturbation 
from the companion star in a binary triggers the same mech- 
anism as that which give rise to the observed distribution 
of eccentricities for planets around single stars. This mech- 
anism is most likely scatterings caused by strong planet- 
planet interactions. For planets which orbit stars that are 
currently single, the strong planet-planet interactions might 
have occurred because the planetary system in which they 
formed was intrinsically unstable (see discussion in section 
1). However, it may be that some of the observed extraso- 
lar planets come from solar-system-like planetary systems. 
In these strong planet-planet interactions occurred because 
the host star suffered a close encounter with another star or 
was exchanged into a binary in an exchange encounter in a 
young stellar cluster. 

Such exchange encounters between single stars and bi- 
nary systems may in fact be rather common in young stellar 
clusters, in which most stars form. We define a singleton 
as a star which did not form in a binary, has never later 
spent time within a binary and has never suffered a close 
encounter with another star. Not all of the stars which are 
single today are singletons. N-body simulations show, that 
as a lower bound, between 5 and 10 per cent of the current 
single field-stars with a mass close to 1 Mq has previously 
either suffered a close enc ounter with another sta r or been 
part of a binary system l|Malmberg et all l2007bl ). From a 
statistical analysis of radial-velocity searches for extrasolar 
planets, it has been estimated that about 7 per cent of all 
solar-mass stars h ave planets on orbits tighter than 5 au 
jMarcv et al.ll2005l l. Hence, external perturbations by other 
stars on solar-system-like planetary systems could account 
for a significant fraction of the observed systems with sepa- 
rations between 1 and 6 au. It is however also possible that 
the contribution is very small, depending on, for example, 
how common solar-system-like planetary systems are and 
how effective fly-bys are for triggering strong planet-planet 
interactions. 

4 SIMULATIONS 

We performed more than 500 simulations of several different 
planetary systems in binaries, using the publicly availabl e 
MERCURY code (|Chambersl fl999l ; [chambers et alj|2002j ). 
All simulations were run for 10 8 years. In all our simula- 
tions we have closely monitored the energy and angular mo- 
mentum conservation, and if it failed using the appropriate 
symplectic integration algorithm (hybrid or wide-binary) in- 
cluded in MERCURY, we re-ran the simulation, using ex- 
actly the same initial conditions, but with a Burlish-Stoer 
algorithm. In the end we used the the Burlish-Stoer method 
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Figure 2. In the upper panel we plot the eccentricity, e, as a 
function of time, t (in years) for four planets inside a binary sys- 
tem. The planets all had a mass equal to that of Jupiter and were 
initially placed on the same orbits as those of the gas giants in the 
solar system. The properties of the companion star were a c = 300 
au, e c = 0.3, i c = 60° and m c = O.GMq. In the lower panel we 
plot the semi-major axis, a (in au) as a function of time for the 
same system. The Kozai Mechanism induces large eccentricities 
in the planetary orbits which leads to strong interactions between 
planets. After only 3 Myr all but one of the planets have been 
ejected, leaving that remaining on a much tighter and much more 
eccentric orbit. 



for the majority of our simulations, since in systems where 
many strong close-encounters between planets occurred, the 
energy conservation using the symplectic algorithms was not 
good enough. 

One woul d expect there to b e a wide variety of plane- 
tary systems dLevison et al.l ll998) but here we divide them 
into two different groups: hierarchical and democratic. We 
only consider solar-system-like planetary systems, i.e. sys- 
tems containing giant planets with separations of around 

5 au or larger. A hierarchical system is dominated by its 
most-massive planet, an example of this is our own Solar 
System. A democratic system on the other hand consists of 
several planets of roughly similar, but not necessarily equal, 
mass. Due to current observational limits, we have not yet 
observed any multiple planet systems with planets outside 

6 au, apart from our own solar system. To simulate a demo- 
cratic system, we used the four giant planets in the solar 
system but set their masses equal to that of Jupiter. To 
check that this system is stable we performed several sim- 
ulations of it around a single star and found no signs of 
any secular trends in the orbital elements of the planets. 
When placed in a binary, we call this system 4 Jupiters-In- 
Binary (4 JIB). Within this system the perturbation of the 
companion star leads to large eccentricities and thus strong 
interactions between the planets, typically resulting in the 
ejection of all but one planet within a few Myr, leaving the 
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remaining planet on a tighter and more eccentric orbit (see 
Fig. 2). 

The effect of the companion star on the planets is to 
slightly perturb the outer planet, leading to strong planet- 
planet interactions in the system. Whether this happens 
or not depends both on the properties of the planetary 
system and on the properties of the companion star. We 
have kept the properties of the companion star constant in 
all the runs from which we generate eccentricity distribu- 
tions below, with a c = 300 au, ra c = O.6M0ande c = 0.3. 
These properties are representative of binaries formed in 
exchange encounters in young stellar clusters (see Fig. 5 in 
iMalmberg et~ai]|2007bl ). Furthermore, the orientation of the 
orbit of the companion star with respect to the orbits of the 
planets in such binaries is random, and thus we assume that 
the inclination of binaries produced in exchange encounters 
is isotropically distributed. 

In order to investigate how the final eccentricity dis- 
tribution is affected by changing the binary properties we 
have also performed simulations of the system consisting of 
four Jupiter-mass planets in a binary, with several differ- 
ent values of the semi-major axis of the binary. These show 
that there are essentially only two outcomes. Either strong 
planet-planet interactions are induced, resulting in an eccen- 
tricity distribution like that of the 4JIB with a c = 300 au 
or the system remains unperturbed, resulting in an eccen- 
tricity distribution like that of the 4 J. The probability that 
strong planet-planet interactions will be induced in the four 
Jupiter system decreases with increasing binary semi-major 
axis. It is close to one for all inclinations up to a c ~ 800 au 
and reaches zero for 1000 < a c < 1500 au, the exact value 
depending on the inclination of the companion star. When 
comparing with the observed eccentricity distribution we 
can thus account for different binary properties by combin- 
ing, for example, the 4J and the 4JIB distribution. Chang- 
ing the mass of the companion star is essentially equivalent 
to changing its semi-major axis, since it only changes the 
strength of the perturbation. Hence we expect the same out- 
come to be true when varying the mass as when varying the 
semi- major axis. 

In order to generate a distribution of eccentricities from 
our simulations we "observe" the eccentricities and semi- 
major axes of all the planets in each simulation at several 
random times. To avoid the initial strong planet-planet in- 
teraction phase, we only considered the second half of our 
simulations, that is, the last 50 Myr. We then discarded all 
the planets with semi-major axes smaller than 1 or larger 
than 6 au (the observational limit) and those with an eccen- 
tricity greater than 0.6. 



5 RESULTS 

We find that the democratic 4JIB system gives rise to a cu- 
mulative eccentricity distribution which goes roughly as e 2 
up to 0.6. We have also performed the same set of simu- 
lations for a system where we instead set the mass of the 
four giants in the solar system equal to that of Uranus. This 
system gives rise to a very similar eccentricity distribution. 
It is interesting to note that the eccentricity distribution 
of democratic systems is very similar to the thermal eccen- 
tricity distribution (/(e) = e 2 ) found for wide stellar bina- 
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Figure 3. The cumulative eccentricity distributions, f(c), gener- 
ated from our simulations of: 4 Jupiter mass planets In a Binary 
(4 JIB); 4 giants of the Solar System In a Binary (SSIB); 4 Jupiter 
mass planets around a single star (4J); a single planet in a binary. 
Also included is the thermal distribution of eccentricities and the 
observed extra— solar planets with separations between 1 and 6 au 
and eccentricity below 0.6. 



rics (|Heggidll975l ). We conclude that democratic systems 
in general produce an excess of eccentric systems compared 
to the observed planets. To compare this result with hier- 
archical systems, we simulated the four giants in the Solar 
System In a Binary (SSIB) and calculated the resulting ec- 
centricity distribution. We find that there is a large excess of 
low-eccentricity systems produced compared to the observed 
planets. 

We have also calculated the evolution of eccentricity 
of a single planet, considering the same stellar binary as 
used earlier, using the semi-anal ytic formulae derived from 
the so-called Kozai Mecha nism 1 Kozail 1962; In nanen et al.l 
1 19971 ; ICarruba et~afl I2002T ) . This system can be thought of 
as an extreme example of a hierarchical system. The result- 
ing cumulative eccentricity distribution is very similar to 
that which we found from our SSIB simulations (see also 
iTakeda fc Rasiol l2005h . Planetary systems containing four 
Jupiters around a single star (4J) give rise to a cumulative 
eccentricity distribution having an excess of low eccentricity 
planets compared to observations. In Fig. 3 we plot the cu- 
mulative distributions of all the above-mentioned systems. 
None of them individually provide a good match with the 
observed distribution. The KS-probability when comparing 
the 4JIB eccentricities with those of the observed planets 
is about 0.001 and when comparing the SSIB eccentricities 
to those of the observed planets the KS-probability is es- 
sentially zero. Hence, the eccentricities of the planets in our 
4JIB and our SSIB samples are very different from the ec- 
centricities of the observed extrasolar planets. 



© 2008 RAS, MNRAS 000, [Hg] 



On the origin of eccentricities among extrasolar planets 5 




0.2 0.4 0.6 0.8 

e 



Figure 4. The cumulative eccentricity distribution, f(e), of the 
observed planets with e < 0.6 compared to that from a combi- 
nation of 30 per cent hierarchical (solar system in binary and 
single planet in binary) and 70 per cent democratic systems (four 
jupitcrs in binary). A KS-test of the two samples gives a proba- 
bility of 15 per cent, which means that they are consistent with 
being drawn from the same population. 



6 DISCUSSION 

Since democratic systems in binaries produce an excess of 
eccentric systems compared to observations and hierarchi- 
cal systems produce an excess of low eccentric systems it 
seems plausible that a combination of the two can provide 
a reasonable match to the observations. We plot an exam- 
ple of such a combination in Fig. 4. In this particular case 
70 per cent of the planets comes from the four jupiters in 
binary (4JIB) systems, while the remaining 30 per cent is 
a mix of solar system in binary (SSIB) and single planet 
in binary systems. Comparing this sample of eccentricities 
with the observed sample using a KS-test gives a probabil- 
ity of 0.15. This is much larger than the probability when 
using a pure sample of either democratic or hierarchical sys- 
tems, and as can be seen in Fig. 4 the match is rather good. 
This fit is surprisingly good considering that we have simply 
taken planetary systems derived from our own solar system. 
This does not suggest that all planetary systems formed are 
either purely democratic or purely hierarchical. However, 
if the contribution to the observed extrasolar planet sam- 
ple from solar-system-like planetary systems exchanged into 
(and sometimes out of) binaries traces the observed eccen- 
tricity distribution, the "average" solar-system-like plane- 
tary systems must be significantly more democratic than 
the four giants in the solar system. Simulations of planetary 
formation shows that a wide variety of planetary systems 
are formed, as can, for example, be seen in the extensive 
catalogue of planeta ry systems produced in simulations by 
iLevison et al.l l|l998t) . It is encouraging to note that a signifi- 
cant fraction of systems formed in such are, by our definition, 
democratic. 



7 SUMMARY 

Most stars form in some sort of cluster or association and 
hence so do most planetary systems. In such crowded places 
initially single stars may be exchanged in and out of binary 
systems and/or pass close to other stars. We have performed 
a large set of simulations of solar-system-like planetary sys- 
tems whose host star have been exchanged into a binary. 
Because the binary was formed in an exchange encounter 
the orientation of the orbits of the planets with respect to 
the orbit of the companion star is completely random. This 
means that the so-called Kozai mechanism operates in many 
of the systems, causing strong planet-planet interactions to 
occur. These in turn lead to the ejection of one or more 
planets, leaving those remaining on tighter and more eccen- 
tric orbits. We find that democratic planetary systems (in 
which the gas giants all have rather similar masses) in bina- 
ries produce an excess of highly eccentric systems compared 
to the observed extrasolar planets while hierarchical plane- 
tary systems (like our own solar system) in binaries produce 
an excess of low-eccentric systems. A combination of hier- 
archical and democratic systems in binaries does however 
provide a good match to the observed eccentricities of ex- 
trasolar planets. 
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